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(§4} Halftoning with enhanced dynamic range and edge enhanced error diffusion. 

(57) Gray pixel values, where the value of the pixel 
has one of a relatively large number of input 
levels i are directed through an error diffusion 
quantization process, for the conversion to one 
of a relatively small number of values. Each 
pixel of the gray level data (18) is modified (10) 
by a correspondingly weighted error correction 
term (20,22,24,28) or terms from previously pro- 
cessed pixels, generating a modified pixel value 
(12). This modified pixel value is compared (14) 
_to a threshold value, the threshold value deter- 
mined aclaptively (15), based on 1) a dither 
pattern, having a dynamic range greater than 1, 
and 2) the value of the pixel to be thresholded, 
and optionally on its predetermined 
neighborhood, multiplied by a value represent- 
ing an edge enhancement factor. The difference 
value between the modified pixel value (12) and 
the output pixel value (16) is distributed in 
accordance with a weighted distribution 
scheme to a set of neighboring unprocessed , 
pixels, increasing or decreasing the gray values '. 
of such pixels. 
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This invention relates to quantizing pixel values in 
an image formed by a plurality of gray level pixels. 

Image information, be it color or black and white, 
is commonly generated in a bitmap format where the 
bitmap comprises a plurality of gray level pixels, i.e. 
pixels that are defined by digital values, each value 
representing a gray level among a number of gray lev- 
els. Thus, in an 8 bit system, 256 levels of gray or 256 
colors are present, where each level represents an in- 
crement of gray between black and white. In the case 
of color bitmaps, where three defining colors or sep- 
arations each incl_ude-256Jevels of information, there 
may be more than 16 million colors defined by a gray 
bitmap. 

Usually, bitmaps in such a gray level format are 
unprintable by standard printers. Standard printers 
print in a limited number of levels, either a spot or a 
no spot in the binary case, or a limited number of lev- 
els associated with the spot, for example, four in the 
quaternary case. Accordingly, it is necessary to re- 
duce the gray level image data to a limited number of 
levels so that it is printed. Besides gray level informa- 
tion derived by scanning, certain processing techni- 
ques such as those described, for example, in our Eu- 
ropean Patent Application No. 91 309 640.0, produce 
gray level pixel values which require conversion to a 
limited set of "legal" or output values. 

One standard method of converting gray level 
pixel values to binary level pixel values is through the 
use of dithering or halftoning processes. In such ar- 
rangements, over a given area having a number of 
gray pixels therein each pixel value of an array of 
gray level pixels within the area is compared to one 
of a set of preselected thresholds (the thresholds are 
stored as a dither matrix and the repetitive pattern 
generated by this matrix is considered a halftone cell) 
as taught, for example, in US-A 4,149,194 to Holla- 
day. The effect of such an arrangement is that, for an 
area where the image is gray, some of the thresholds 
within the dither matrix wilt be exceeded, i.e. the im- 
age value at that specific location is larger than the 
value stored in the dither matrix for that same loca- 
tion, while others are not. In the binary case, the pix- 
els or cell elements for which the thresholds are ex- 
ceeded might be printed as black, while the remaining 
elements are allowed to remain white, dependent on 
the actual physical quantity described by the data. 
The effect of the distribution of black and white over 
the halftone cell is integrated by the human eye as 
gray. Dithering or halftoning presents problems, how- 
ever, in that the amount of gray within an original im- 
age is not maintained exactly over an area, because 
the finite number of elements inside each dither ma- 
trix - and therefore halftone cell - only allows the re- 
production of a finite number of gray levels, i.e. equal 
or less than the number of elements in the cell plus 
one, or less. The error arising from the difference be- 
tween the output pixel value and the actual gray level 



pixel value at any particular cell is simply thrown 
away. This results in a loss of image information. In 
particular, dithering introduces coarse quantization 
artifacts which are visible in the image areas where 

5 the scene has little variation. This is also known as 
"banding", and is caused by the limited number of out- 
put gray levels available. The "banding" artifacts gen- 
erally increase with decreasing cell size, which is 
identical to a decrease in the number of levels that 

10 can be represented by the halftone cell. 

Algorithms that convert gray images to binary or 
other number of level images whife^attempting to pre- 
serve local density exist, and include among them er- 
ror diffusion, as taught, for example, in "An Adaptive 

15 Algorithm for Spatial Greyscale" by Floyd and Stein- 
berg, Proceedings of the SID 17/2, 75-77 (1976) 
(hereinafter, "Floyd and Steinberg"). Another, more 
elaborate method would be the error diffusion tech- 
nique of US-A 5,045,952 to Eschbach, which serves 

20 to provide image dependent edge enhancement, as- 
signed to the same assignee as the present inven- 
tion. Additional modifications to the error diffusion al- 
gorithm as taught by Floyd and Steinberg have been 
proposed, e.g.: a different weighting matrix, as 

25 taught, for example, in "A Survey of Techniques for 
the Display of Continuous Tone Pictures on Bilevel 
Displays" by Jarvis et al.. Computer Graphics and Im- 
age Processing, Vol. 5., pp. 13-40 (1976), and in 
"MECCA - A Multiple-Error Correction Computation 

30 Algorithm for Bi-Level Image Hardcopy Reproduc- 
tion" by Stucki, IBM Res. Rep. RZ1060 (1981). Mod- 
ifications of the error calculation and weight alloca- 
tion have been taught, for example, in US-A Patent 
4,924,322 to Kurosawa et. al., US-A 4,339,774 to 

35 Temple, and US-A4,955,055 1 to Ulichney. Other error 
diffusion methods include those described in our Eu- 
ropean Patent Applications Nos. 91 309 640.0 and 92 
308 074.1. 

Error diffusion attempts to maintain gray by mak- 

40 ing the conversion from gray pixels to binary or other 
level pixels on a pixel-by : pixel basis. The procedure 
examines each pixel with respect to a threshold or set 
of thresholds, and the difference between the gray 
level pixel value and the output value is forwarded to 

45 a selected group or set of neighboring pixels, in accor- 
dance with a weighting scheme. 

Another modification to the method error diffu- 
sion was taught by Billotel-Hoffmann and Bryngdahl, 
Proceedings of the SID, Vol. 24/3, (1983), pp. 253- 

so 258 (hereinafter, BrlloteMHoffmann and Bryngdahl). A 
dither matrix is used as a threshold for error diffusion, 
to alleviate the problems of undesired patterns gen- 
erally produced by the error diffusion algorithm. Not- 
ably, however, in Billotet-Hoffmann and Bryngdahl, 

55 the maximum dynamic range of the proposed thresh- 
old modification is set to a. = 1 (page 257 in the afore- 
mentioned reference). 

A problem noted with the use of the standard er- 
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ror diffusion algorithms for printing applications is the 
production of large numbers of isolated black and/or 1 
white pixels which are non-printable by many types of 
printers. The algorithm taught by Bil I otet- Hoffmann 
and Bryngdahf does not improve the printability of the 
general error diffusion algorithm, A method to over- 
come the printability problem is taught by US-A 
4,654,721 to Goertzel t where a method is shown to 
convert a continuous tone image to a bilevel pixel im- 
age. The total error generated in one halftone cell is 
distributed to a predetermined number of adjacent 
halftone cells. !n this way, printable images are gen- 
erated, while the banding artifact is reduced, by alter- 
nating between fixed output dot patterns. Because of 
an inherent lack of partial dots in this process, evi- 
denced as a loss in sharpness, edge detection and 
sharpening was included. See, also, "Digital Halfton- 
ing in the IBM 4250 Printer" by Goertzel etal, (Goert- 
zel) , IBM J. Res. Develop,. Vol 31, No. 1 f January, 
1987. 

It is an object of the present invention to provide 
an improved method and apparatus for quantizing 
pixel values. 

According to one aspect of the invention, there is 
provided a method of quantizing pixel values in an 
original image formed by a plurality of pixels, each 
pixel representing an optical density of the image at 
a location within the image, arid having an original 
optical density value associated therewith selected 
from one of a set of c original optical density values 
that has a number of members larger than a desired 
output set of d optical density values, and having ady- 
namic range associated therewith, the steps com- 
prising for each pixel, adding a previously determined 
error term to the original optical density value thereof, 
to derive a modified pixel value; determining for the 
pixel to be quantized, a threshold level that is a func- 
tion of the original optical density of the pixel, amd op- 
tionally its predetermined neighborhood, to be quan- 
tized, and a threshold value selected from a set of 
threshold values having a dynamic range that is 
greater than that of the original image; applying the 
determined threshold level to each optical density 
value of each pixet in the image, to produce a output 
optical density value that is a member of the desired 
output set, each member of the output set of d optical 
density values a legal output value; determining an 
error term that is a difference between the output opt- 
ical density value and the modified optical density val- 
ue of the pixel; applying a proportional amount of the 
determined error term to the original optical density 
value of each of a predetermined set of neighboring 
pixels. 

According to another aspect, the invention pro- 
vides an arrangement for quantizing pixel values in an 
image formed by a plurality of pixels, each pixel rep- 
resenting an optical density of the image at a location 
within the image, and having an original optical den- 



sity value associated therewith selected from one of 
a set of c original optical density values that has a 
number of members larger than a desired output set 
of d optical density values, and having a dynamic 

5 range associated therewith defined as 1 , the quantiz- 
ing arrangement comprising means for inputting an 
image pixel to be quantized having an original optical 
density value; means for adding an error term derived 
from the quantization of at least one previous pixel to 

10 the original optical density value to produce a modi- 
fied image pixel; means for determining for the pixel 
to be quantized, a threshold level, that is a function 
of the original optical density of the pixel, and option- 
ally its predetermined neighborhood, to be quantized 

15 and a threshold value that is selected from a set of 
threshold values having a dynamic range that is 
greater than that of the original image; means for 
thresholding the optical density value of each modi- 
fied pixel in the image with the threshold level deter- 

20 mined at said determining means, to produce a new 
optical density value that is a member of the desired 
output value set each member of the output set of d 
new optical density values a legal output value which 
may be represented by an output device; means for 

25 determining an error term that is a difference be- 
tween the new optical density value and the modified 
optical density value; means for storing a set of 
weighted error terms to be applied to each original 
optical density value in a predetermined set of neigh- 

30 boring pixels, each weighted errorterm a preselected 
portion of the error term. 

The invention provides a quantizing method 
wherein a pixel def ined by one of a large number of 
possible levels of optical density is redefined at one 

35 of a smaller number of levels of optical density, using 
an edge enhancing error diffusion algorithm with a 
threshold set in accordance with a dither matrix of 
large dynamic range. 

In one preferred form of the invention, gray pixel 

40 values, where the value of the pixel has one of a rel- 
atively large number of input levels, are directed 
through an error diffusion quantization process, for 
the conversion to one of a relatively small number of 
values. Each pixel of the gray level data is modified 

45 by a correspondingly weighted error correction term 
or terms from previously processed pixels, generat- 
ing a modified pixel value. This modified pixel value 
is compared to a threshold value, the threshold value 
determined adaptively, based on 1) a value given by 

so the dither matrix, having a dynamic range greater 
than 1 (with the possible dynamic range of the input 
image being defined as 1 ), and 2) the value of the pix- 
el to be thresholded, multiplied by a constant value 
representing an edge enhancement factor or by any 

55 of the edge enhancement schemes for error diffu- 
sion. The difference value between the modified pixel 
value and the output pixel value is considered the er- 
ror and is distributed in accordance with a weighted 
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distribution scheme to a set of neighboring unpro- 
cessed pixels, increasing or decreasing the gray val- 
ues of such pixels. The term 'optical density' will be 
used hereinafter to describe the set of values com- 
prising the image data, be it density, intensity, light- 
ness or any other physical quantity described by im- 
age data representation. 

In accordance with another aspect of the inven- 
tion, there is provided a method of quantizing pixel 
values in an image formed by a plurality of pixels, 
each pixel representing an optical density of the im- 
age at a location wit hi n th e image, and having an orig- 
inal optical density value associated therewith select- 
ed from ore of a set of c original optical density values 
that has a number of members larger than a desired 
output set of d optical density values, the steps com- 
prising: adding a sum of error terms derived from the 
quantization of at feast one previous pixel and the 
original optical density value to derive a modified den- 
sity value reflecting the added error applying a 
threshold decision to" each modified optical density 
value of each pixel in the modified original image, to 
produce a new optical density value that is a member 
of the desired output set, each member of the output 
set of d optical densityvaluesalegal output value, the 
threshold decision derived as a function of the input 
image, and a set of dither matrix threshold values 
having a dynamic range greater than the dynamic 
range of original optical density values; determining 
an error term that is a difference between the new 
output optical density value and the modified optical 
density value; applying a proportional amount of the 
error term to the original optical density value of each 
of a predetermined set of neighboring pixels, taking 
into account the fractional error allocations for said 
pixels generated at other pixels than the current one. 

The present invention teaches a method thatcaf- 
cutates error diffusion on a pixel- by-pixel basis, but si- 
multaneously achieves printability by introducing a 
dither matrix into the process, but with a dynamic 
range clearly exceeding the dynamic range of the in- 
put, in contrast to the method taught by Billotet- Hoff- 
mann and Bryngdahl. Additionally, use is made of the 
method described in U3-A 5,045,952 to Eschbach, to 
maintain image sharpness. There is a clear prefer- 
ence for a pixel-by-pixel processing enabled by this 
method, because it allows a better representation of 
image detail on the pixel level. 

These and other aspects of the invention will be- 
come apparent from the following descriptions to il- 
lustrate a preferred embodiment of the invention read 
in conjunction with the accompanying drawings in 
which: 

Figure 1 is a block diagram of a basic system for 
carrying out the present invention; 
Figure 2A shows an example standard dither ma- 
trix (dynamic range = 1), while Figure 2B shows 
the dither matrix modified to have an increased 



dynamic range, here 4; 

Figures 3A and 3B show a flow chart of the de- 
scribed embodiment; 

Figure 4 is a different depiction of the basic sys- 

5 tern described in Figure 1; and 

Figures 5, 6, 7 A, 7B, 8A and SB are an illustrative 
example of the difference between Floyd and 
Steinberg, Billotet- Hoffmann and Bryngdahl, and 
the present invention. 

10 Referring now to the drawings where the show- 

ings areforthe purpose of describing an embodiment 
of the invention and not for limiting same, a basic sys- 
tem for carrying out the present invention is shown in 
Figure 1. In the present case, gray level image data 

15 from image input 8 may be characterized as image 
data, each pixel of which is defined at a single level 
or optical density in a set of optical density levels, the 
number of members in the set of levels being larger 
than desired. Each pixel will be processed in the man- 

20 ner described hereinbelow, to redefine each pixel in 
terms of a new, smaller set of levels. Here, color data 
may be represented by a number of independent 
channels which are handled independently, or the 
color data might be represented as vector data in a 

25 predefined color space, e.g.; RGB, CIELab etc., being 
submitted to vector operations in the thresholding, er- 
ror calculation and correction. One common case of 
this method includes the conversion of data from a 
relatively large set of gray level values to one of two 

30 legal or allowed bin values for printing in a binary prin- 
ter. Another case of this is the conversion of data from 
a relatively large set of color data expressed as red, 
green and blue, or cyan, magenta, yellow and black, 
to one of five legal bin values for printing, as descri- 

35 bed in our European Patent Application No. 92 308 
074,1. 

An input image of the type to be processed as' 
hereinafter described may be represented by a set of 
gray values (gray level pixels) arranged in an array of 

40 £ lines, each line containing n gray values with depth 
b. Gray values are typically expressed as integers, 
with one example falling in the range from 0 to 255, 
although greater or lesser number of levels, as well as 
non-integer representations, are possible. An output 

45 image is considered to consist of pixels, each pixel 
corresponding to an output element that is printed by 
a digital printer or display. A dither matrix, as used 
herein, refers to a set of thresholds defined at loca- 
tions corresponding to pixels over a given area of the 

50 image. Adither matrix is generally smaller than the to- 
tal image and will be replicated in a predetermined 
scheme in order to cover an area of the image. A 
method for an efficient representation of variable an- 
gle halftone cells by a dither matrix and a correspond- 

55 ing replication scheme, is given in US-A 4,149,1 94 to 
Holladay. The output of a process using a dither ma- 
trix is a set of pixel values, having a number of mem- 
bers less than the input set of values. Commonly, the 
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set of output values is binary, either biack or white, or 
a spot or no spot, although the values might be gray. 
The output of the dithering, or the set of pixels that 
are either black or white, together form a "dot". Single 
pixels, black or white, surrounded respectively by 
white or black pixels, are difficult to print with electro- 
photographic devices. For this reason, standard dith- 
er matrices for electrophotographic applications tend 
to cluster pixels together, with a growth pattern that 
begins in a centra! area of the halftone cell and grows 
as more elements of the cell are black. Such a dot pat- 
tern is printable on electrophotographic devices. Dith- 
er matrices are sometimes also referred to as half- 
tone screens. 

With reference to Figure 1 , a stored array of input 
image data or pixels 6, which may be any source of 
image data, directs input image I into the system on 
a pixel- by-pixel basis, where n, £ represents the pos- 
ition of a single pixel l n ,r in the stream of image data. 
ln T f refers in this description to both the pixel that is 
positioned, at nj in the image stream, and the inten- 
sity of the pixel at position n,£. Each input pixel has a 
corresponding error term or value e added to the input 
value t„,f at adder 10, where e n ,f is a sum of error val- 
ues of previous pixels to be directed to l n , fr resulting 
in a modified image, represented by modified pixel 
values, temporarily stored at block 12. The modified 
image, the sum of the input value and the error value 
of previous pixels (l n ,f + £ n , f ), is passed to threshold 
comparator 14- The modified image is compared to 
threshold value(s) T^f to determine an appropriate 
output value B ntf for pixel l„ r such as, for examples, 
in a binary output image, a spot or no spot. The ap- 
plication of variable thresholds is represented by the 
threshold source 15, which will be further explained 
herelnbelow. Once output pixel B n ^ is determined and 
directed to output image storage 16 for eventual trans- 
fer to an output along line 18, the value of B n ,f is sub- 
tracted from the modified image value (l n ,f+En r f) to 
generate an error level £ m from pixel l n#f . The subtrac- 
tion operation is represented by the change sign 
block 20 and subsequent adder 22, with s m represent- 
ing the difference between the modified image value 
(l n ,f +£n r r) and the output value B n j stored to errorfrac- 
tion determination block 24, where weighted portions 
of error e m are calculated, and will be used for updat- 
ing the error delay buffer 28. 

Determination of the threshold at threshold 
source 15 is made in accordance with the method of 
edge enhanced error diffusion taught in US-A 
5,045,952 to Eschbach. That description provides a 
variable threshold level T nt f determined by first calcu- 
lating a modifier based on the input value \„ r £ of each 
pixel and/or its predetermined neighborhood as de- 
scribed in US-A 5,045,952, as represented in the 
threshold modifier block 40. For the simple case of 
determining the modifier based on the input pixel t n f 
without a neighborhood, the modifier value l n ,f x (K- 



1) is then subtracted from the nominal threshold val- 
ue M n r, at threshold calculation 42, to determine 
threshold level T^p to be applied to threshold compar- 
ator block 20, with the e nhan cement factor K to be se- 

5 lected as hereinafter described. Alternatively, it is 
possible to accomplish an equivalent alteration of the 
threshold through the addition of the modifier value 
hrtK" 1 ) to the modified image value (In+E^.,) while 
maintaining the threshold value at level lA^j. 

10 In accordance with the invention, threshold value 

M n f is determined in accordance with the position of 
pixel l nr f in the data stream, represented in Figure 1 
by the clocking of pixel and scan iines commonly gen- 
erated in image processing systems, and a halftone 

15 screen or dither matrix M of threshold values. A set of 
values M stored in a dither matrix storage or memory 
44 is directed on a pixel by pixel, and scan line by scan 
line basis to threshold calculation block 42, through 
the dynamic range multiplier 46, resulting in the 

20 threshold T^f being a superposition of dither matrix 
and input image information. A changing pattern of 
threshold values is therefore superimposed over the 
edge enhancing error diffusion algorithm. 

In accordance with another aspect of the inven- 

25 tion, the threshold values of dither pattern M are se- 
lected having a dynamic range greater than the pos- 
sible dynamic range of the input pixel values. Using 
a standard dither pattern with a dynamic range that 
is equal to or less than the range of the input values, 

30 such as perhaps that shown in Figure 2A for the case 
of input values in the range from 0 to 255, each 
threshold level in the dither pattern is multiplied by a 
factor /, where / is much greater than 1 and increas- 
es with halftone cell or dot size. For example, for a 4 

35 x 2 screen, / = 4; for a 8 x 4 screen, / = 8; and for a 
10x5 screen, / = 10. These values are only exam- 
ples, and even for the particular screen sizes cited, 
other values of / may be provided. It will be under- 
stood that the value of / is greater than 1, and in- 

40 creases with dot size. Of course, screens with such 
threshold values could be constructed, butuserstend 
to have a large selection of screens in their posses- 
sion. Accordingly, a simple multiplier arrangement, 
such as that provided by the combination of dither 

45 matrix storage 44 and dynamic range multiplier 46, is 
desirabie. As previously noted, this value /M n , f is 
then used in the calculation of T ni f. Alternatively, the 
original dfther matrix (Figure 2A) may be stored into 
the dither storage block 44 through a multiplier. In 

so such a case, block 46 can be eliminated and the ex- 
tended dynamic range dither matrix (Figure 2B) can 
directly be addressed by the pixel and scan line dock. 
In any case, the value directed to threshold calculator 
operates as if the example dither matrix of Figure 2B 

55 were being used. 

Enhancement factor K may be held constant or 
may vary as a function of the input image content, lo- 
cal or global, within the continuous tone input image. 
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In this particular application, K is selected to be about 
trie same value as the dynamic range multiplier used 
for the dither matrix, with another embodiment using 
a K larger than the dynamic range multiplier of the 
dither matrix. As will become apparent, K is a factor 
in determining how closely image detail is preserved. 
Accordingly, a larger K tends to improve the apparent 
sharpness of the image. 

With reference to Figure 4, there is provided a dif- 
ferent representation of the present invention, where 
the error compensation is not done via a modified im- 
agers shown in Figure 1 , but by using a continuously 
updated threshold adjustment delay buffer. Accord- 
ingly, a stored array of input image data or pixels 100, 
which may be any source of image data, directs input 
image I into the system on a pixel-by-pixel basis. The 
image is passed to threshold comparator 114. The im- 
age is compared to threshold value(s) T n j to deter- 
mine an appropriate output value B nJ for pixel l n , ? t 
such as, for example, in a binary output image, a spot 
or no spot, which is stored to output buffer 116, The 
output value B n(P is changed in sign, and added at ad- 
der 122 to the original value of pixel l n , P , so that a 
threshold adjustment allocation 0^=^ -B^. Thresh- 
old adjustment allocation D m represents an error in 
the thresholding process, and a weighted distribution 
thereof to a selected group of neighboring; unpro- 
cessed pixels. The original value of pixel plus op- 
tionally its neighborhood, is also used forthe creation 
of a modifier value based on the input value \ n j of 
each pixel, as represented in the threshold modifier 
block 140. A modifier value l n ? x(K-1) is then subtract- 
ed from the expanded dynamic range threshold value 
/M af , at threshold calculation 142, to determine 
threshold level T n (> to be applied to threshold compar- 
ator block. The dither matrix values /M n ^ are derived 
as previously described, and represented as extend- 
ed dynamic range dither matrix storage 144. At 
threshold adjustment delay buffer 128, threshold lev- 
el T n ,f is added to the sum of threshold adjustment al- 
locations D n( p directed thereto, for the creation of a 
new threshold T 0f f. 

With reference now to Figure 3, there is shown 
a flow chart demonstrating the steps of the inventive 
process of quantizing pixel values in an image formed 
by a plurality of pixels, each pixel representing an opt- 
ical density of the Image at a location within the im- 
age, and having an original optical density value as- 
sociated therewith selected from one of a set of c orig- 
inal optical density values that has a number of mem- 
bers larger than a desired output set of d optical den- 
sity values, where steps 110, 120, 130, 140, 150, 160, 
170 and 180 area standard error diffusion method of; 
110) receiving an input pixel having a value i n( p; 120) 
adding an error term e nt ? from the error delay buffer 
to the pixel value \ n j\ 130) storing the value l Ri f + e n> f 
as a modified pixel value; 140) thresholding the modi- 
fied pixel value with T P| P to obtain a new output value 



B n ^ and directing B n j to an output; 150) changing the 
sign of output value B nfP 160) adding -B nif to l n< r + £ n ,r 
to obtain quantizing error e^, 170) determining the 
weighted error for each of the currently unprocessed 

5 pixels to which error in the quantization of pixel l n){ ' is 
to be directed, and 180) updating an error storage or 
buffer holding the errors that are added to the future 
pixels at step 120. Here, the unprocessed pixel set 
r f • 'n- i,p*i, h.f+i, In+i ,f*i} shown in block 170 and 180 

10 is a common set for error diffusion, but other sets are 
also possible and are intended to be included in this 
description. In accordance with the invention, and 
shown in Figure 3B, at step 140, an additional thresh- 
old determination is provided to determine threshold 

15 level T rtrf , including onsubstep 142, receiving l n?P ; sub- 
step 144, calculating a modifier based on the input 
value l ni f of each pixel and an enhancement factor K; 
substep 146, subtracting from threshold value M rt|f , 
the threshold modifier value 1 rtf (K-1) to determine 

20 threshold value T nt r to be applied to threshold the 
thresholding step 140, It can also be appreciated that 
alternatively subblock 142 receives the input optical 
density values of a set of pixels including l nP and that 
the threshold modifier in subblock 144 is calculated 

25 dependent on this set, A simple example is to use a 
filtered version of the input image as input to subblock 
144 as taught in US-A 5,045,952. At step 148, assum- 
ing a standard dither matrix, each value of the dither 
matrix M is multiplied by a value / which is larger than 

30 1 , and increases in size with the number of elements 
of the dither matrix, so that the output of substep 
148 reflects the desired expanded dynamic range of 
the dither matrix values. 

With reference now to Figure 5, there is provided 

35 an example of the inventive method. An array of input 
optical density values is shown in Figure 6, where 
white = 255 and black = 0. The possible input dynamic 
range is assumed to be 0 to 255, i.e. 256 levels. For 
this example, consider the dynamic range to be de- 

40 fined as being 1 . The image data represents a change 
over from an area of darker gray (pixel value = 70) to 
lighter gray (pixel value = 200) being a segment of a 
larger image. Figure 6 shows the output optical den- 
sity values generated by the method described by 

45 Floyd and Steinberg. The large number of isolated 
black or white pixels is apparent. See, for example, 
pixel #4 in scanline #3, where a value of 255 is sur- 
rounded by 0's, or pixel #13 in scanline #12, where a 
value of 0 is surrounded by 255's. Figure 7B shows 

50 the method described by Billotet- Hoffmann and 
Bryngdaht, using a 32 element dither matrix shown in 
Figure 7A (describing a 33 level halftone cell at 45 de- 
grees inclination) with a maximum dynamic range 
factor of "1 The output optical density values are dif- 

55 ferent from Figure 6, but a large number of isolated 
pixels still remain. Figure 8B shows the output optical 
density values of the inventive method, using, for this 
example, the 32 element dither matrix shown in Fig- 
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ure 8A, a dynamic range factor of * r 8" and a threshold 
modifier factor of K= 12, It is apparent from Figure SB 
that the image output structure has fewer individual 
black and white pixels. . 

The invention has been described with reference 
to a particular embodiment. Modifications and alter- 
ations will occur to others upon reading and under- 
standing this specification. It should be noted that the 
use of the term "thresholding" throughout this inven- 
tion is meant to encompass other ways of making a 
selection of the output optical density value. It should 
be noted that the present invention can be equiva- 
lent ly described by the characterization of Figure 4 
and it is intended that all such modifications, altera- 
tions and changes in notation are included insofar as 
they come within the scope of the appended claims 
or equivalents thereof. 



cell. 

3. The method as defined in claim 1, wherein the 
function determining the threshold level is the 

5 difference between a product of the original opt- 

ical density value of the pixel to be quantized and 
a selected enhancement factor K -1, and a 
threshold value selected from the set of threshold 
values having a dynamic range that is greater 
10 than that of the original image. 

4. The method as defined in claim 3, wherein the 
enhancement factor K is selected to be greater 
than or equal to the dynamic range of the set of 

15 threshold values. 

5. The method as defined in claim 1 , where the dy- 
namic range of the set of threshold values is be- 
tween 4 to 8 times the dynamic range of the orig- 
inal image. 

6. An arrangement for quantizing pixel values in an 
image formed by a plurality of pixels, each pixel 
representing an optical density of the image at a 
location within the image, and having an original 
optical density value associated therewith select- 
ed from one of a set of c original optical density 
values that has a number of members larger than 
a desired output set of d optical density values, 
and having a dynamic range associated there- 
with defined as 1, the quantizing arrangement 
comprising: 

means for inputting an image pixel to be 
quantized having an original optical density val- 
ue; 

means for adding an error term derived 
from the quantization of at least one previous pix- 
el to the original optical density value to produce 
a modified image pixel; 

means for determining for the pixel to be 
quantized, a threshold level, that is a function of 
the original optical density of the pixel, and op- 
tionally its predetermined neighborhood, to be 
quantized and a threshold value that is selected 
from a set of threshold values having a dynamic 
range that is greater than that of the original im- 
age; 

means for thresholding the optical density 
value of each modified pixel in the image with the 
threshold level determined at said determining 
means, to produce a new optical density value 
that is a member of the desired output value set, 
each member of the output set of d new optical 
density values a legal output value which may be 
represented by an output device; 

means for determining an error term that 
is a difference between the new optical density 
value and the modified optical density value; 



Claims 20 

1. A method of quantizing pixel values in an original 
image formed by a plurality of pixels, each pixel 
representing an optical density of the image at a 
location within the image, and having an original 25 
optical density value associated therewith select- 
ed from one of a set of c original optical density 
values that has a number of members larger than 
a desired output set of d optical density values, 
and having a dynamic range associated there- 30 
with, the steps comprising: 

for each pixel, adding a previously deter- 
mined error term to the original optical density 
value thereof, to derive a modified pixel value; 

determining for the pixel to be quantized, 35 
a threshold level that is a function of the original 
optical density of the pixel, and optionally its pre- 
determined neighborhood, to be quantized, and 
a threshold value selected from a set of threshold 
values having a dynamic range that is greater 40 
than that of the original image; 

applying the determined threshold level to 
each optica! density value of each pixel in the im- 
age, to produce a output optical density value 
that is a member of the desired output set, each 45 
member of the output set of d optical density val- 
ues a legal output value; 

determining an error term that is a differ- 
ence between the output optical density value 
and the modified optical density value of the pix- so 
el; 

applying a proportional amount of the de- 
termined error term to the original optical density 
value of each of a predetermined set of neighbor- 
ing pixels. 55 

2. The method as defined in claim 1, wherein the 
the set of threshold values describes a halftone 
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means for storing a set of weighted error 
terms to be applied to each original optical den- 
sity value in a predetermined set of neighboring 
pixels, each weighted error term a preselected 
portion of the error term. 

7. The arrangement as defined in claim 6, wherein 
the the set of threshold values describes a half- 
tone cell. 



10 



8. The arrangement as defined in claim 6, wherein 
the function determining the threshold level is a 
difference between a product of the original opt- 
ical density value of the pixel to be quantized and 
a selected enhancement factor K-1, and the 15 
threshold value selected from a set of threshold 
values having a dynamic range that is greater 
than that of the original image. 

9- The arrangement as defined in claim 8, wherein 20 
the enhancement factor K is selected to be great- 
er than or equal to the dynamic range of the set 
of threshold values. 

10, The arrangement as defined in claim 6, where 25 
the dynamic range of the set of threshold values 
is between 4 to 8 times the dynamic range of the 
original image. 
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(54) Halftoning with enhanced dynamic range and edge enhanced error diffusion. 

(57) . Gray pixel values, where the value of the pixel 
has one of a relatively large number of input 
levels, are directed through an error diffusion 
quantization process, for the conversion to one 
of a relatively small number of values. Each 
pixel of the gray level data (18) is modified (10) 
by a correspondingly weighted error correction 
term (20,22,24,28) or terms from previously pro- 
cessed pixels, generating a modified pixet value 
(12). This modified pixel value is compared (14) 
to a threshold value, the threshold value deter- 
mined adaptively (15), based on 1) a dither 
pattern, having a dynamic range greater than 1, 
and 2) the value of the pixel to be thresholded, 
and optionally on its predetermined 
neighborhood, multiplied by a value represent- 
ing an edge enhancement factor. The difference 
value between the modified pixel value (12) and 
the output pixel value (16) is distributed in 
accordance with a weighted distribution 
scheme to a set of neighboring unprocessed 
pixels, increasing or decreasing the gray values 
of such pixels. 
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